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Abstract.  In  this  paper,  the  problems  of  constructing  a  mathematical  model  for  small 
linear  oscillations  of  an  orthotropic  plate  carrying  concentrated  masses  in  a  finite  set 
of  points  connected  with  a  fixed  base  by  elastic  struts  (springs)  with  certain  stiffness 
coefficients,  as  well  as  the  direct  and  inverse  spectral  problems  for  such  a  model  in  a  new 
formulation  are  considered. 

Key  Words  and  Phrases:  mathematical  model,  oscillation,  orthotropic  plate,  direct 
spectral  problem,  inverse  spectral  problem. 

2010  Mathematics  Subject  Classifications:  35P05,  47F05 


1.  Introduction 

In  this  paper,  the  mathematical  models  describing  vibrations  of  elastic  plates 
and  shells  are  studied.  Plates  are  equipped  with  additional  fasteners  and  con¬ 
tain  concentrated  masses  at  some  points.  At  the  same  time,  structural  elements 
(masses,  fastenings)  have  exact  and  relatively  small  dimensions  compared  with 
the  plate.  Similar  designs  are  common  in  various  fields  of  mechanics  and  engi¬ 
neering. 

The  frequencies  and  modes  of  oscillations  significantly  depend  on  the  elas¬ 
ticity  coefficients  and  point  masses  of  fasteners,  which  entails  the  possibility  of 
controlling  these  characteristics  using  the  available  fastener  parameters.  In  ad¬ 
dition  to  the  control  problems  for  such  systems,  the  problems  of  diagnosing  (or 
identifying)  the  parameters  of  the  connections  from  known  data  (the  finite  set  of 
eigenvalues  of  oscillations)  are  relevant. 

In  this  paper,  a  different  approach  to  the  study  of  such  structures,  based 
on  the  theory  of  multiparameter  inverse  spectral  problems  (briefly  MPISP)  are 
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proposed.  For  the  MPISP,  it  is  appropriate  to  apply  the  methods  and  approaches 
developed  in  [1,  2,  8],  which  differ  significantly  from  the  methods  of  classical 
inverse  spectral  problems. 

Let  us  directly  proceed  to  the  mathematical  model  of  the  structure.  Let 
an  orthotropic  plate  fl  be  given  carrying  the  concentrated  masses  rn3  at  certain 
points  Mj  =  Mj  (xj ,  y3 ) ,  which  are  connected  to  the  fixed  base  by  elastic  stiffness 
springs  kj.  We  will  assume  the  validity  of  the  Kirchhoff-Love  hypotheses  with 
the  boundary  conditions  of  any  of  three  possible  types  (free  edge,  hinge  bearing, 
pinching).  The  interaction  of  point  masses  with  elastic  columns  (springs)  is  well 
managed  to  be  modeled  using  Dirac  delta  functions. 

It  is  known  that  within  the  assumptions  on  the  model,  the  natural  oscillations 
of  the  plate  are  described  by  a  biharmonic  operator.  In  our  case,  this  operator  has 
distribution  coefficients  of  5  -like  form.  Such  operators  were  studied,  for  example, 
in  [7,  3], 

Below  it  will  be  shown  that  the  form  of  the  plate’s  natural  oscillations  is 
described  by  an  operator  pencil.  A  direct  problem  is  to  investigate  the  nature 
of  the  spectrum  of  this  operator  pencil.  The  inverse  spectral  problem  consists  in 
determining  the  coefficients  of  spring  stiffness  and  the  values  of  the  concentrated 
masses  from  the  known  finite  set  of  natural  frequencies  of  oscillations. 


2.  Construction  of  a  mathematical  model 

We  assume  that  at  rest  the  form  of  the  considered  plate  is  a  simply  connected 
bounded  domain  fl  C  R2  with  a  smooth  boundary  L  =  dfL  Mj  =  ( Xj,yj )  € 
are  the  points  where  elastic  supports  with  elasticity  coefficients  kj  and  point 
masses  rrij ,  j  =  1 ,  n  are  found. 

Let  the  function  W(M,t),  where  M  =  M(x,y)  €E  D,  be  the  magnitude  of  the 
deviation  from  the  equilibrium  state,  and  either  the  free  edge  condition,  or  the 
hinged  support  condition,  or  the  pinch  condition  be  satisfied  on  the  boundary. 
Without  loss  of  generality,  we  shall  further  assume  that  for  W (M,  t )  the  condi¬ 
tions  of  ’’pinching”  are  fulfilled:  W |r  =  0,  yj)p|r  =  0.  Denote  these  conditions  by 
l(W)  |r  =  0. 

The  plate  thickness  is  assumed  to  be  h ,  h  «  1,  density  is  denoted  by  p , 
Young’s  modulus  by  E  and  Poisson’s  ratio  by  a  .  The  potential  energy  of  the 
plate,  together  with  the  energy  of  the  elastic  supports,  has  the  following  form: 


n  =  D  ((Wxx  +  Wyyf  -  2(1  -  a)(WXXWyy  ~  W^))  dXdV  + 
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1  ” 

+  9  ^2kjW(xj,yj)2, 


(1) 


3= 1 


where  D  =  a  cylindrical  rigidity  of  the  plate.  The  kinetic  energy  of 

the  entire  plate,  taking  into  account  the  added  masses,  has  the  following  form: 


^  [  \Wf\dxdy  +  ^^mj\Wt{xj,yj,t)^ 


T=  — 

In  this  case,  the  action  functional  has  the  form 

t-2 


S=  (T-U)dr 


(2) 


(3) 


Considering  S  on  the  set  of  harmonic  oscillations  of  the  form 

W(x,y,t)  =  elLVtu(x,y), 
we  obtain  a  quadratic  functional 

S(u,u)  =  D  f  (|  Au|2  -  2(1  -  <j)(uxxuyy  -  ulyj)  dxdy— 


a 


j  \u2\dxdy  +  ^  ^2(kj  -  u2mj)\u(xj,  yj)\2. 
o  i=l 


(4) 


Note  that  the  problem  of  finding  the  natural  oscillations  of  a  system  is  equiv¬ 
alent  to  the  search  for  the  extremals  of  quadratic  form  S(u,  u),  see  e.g.  [6]. 

Now,  consider  the  problem  in  dimensionless  quantities,  and  also  introduce  the 
divergent-gradient  operator 


Aq  u  =  div(a\/)u, 


where 


a  = 


an  ai2 
«2i  a  22 


a*  =  a  >  0. 

Let’s  pass  from  (4)  to  the  quadratic  form  in  dimensionless  quantities 
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(5) 

As  the  domain  of  the  quadratic  form  Q{u,u),  it  is  natural  to  consider  the  space 

DQ  =  {ue  W2(n)\l(u)\dn=o}. 

Let  us  introduce  the  Hilbert  space  H\  with  the  norm 

I!  u  llffi=  /( |A0u|2  +  \u\2)dxdy. 

n 

Then  it  is  obvious  that  the  norms  ||  •  \\h1  and  |j  •  \\W2  are  equivalent. 

Thus,  the  problem  of  finding  the  natural  vibrations  of  a  plate  is  reduced  to 
finding  extremals  (eigenfunctions)  of  the  quadratic  functional  Q(u,u,  A)  in  the 
class  of  functions  Dq. 

Let  Q(u,v,  A)  be  the  sesquilinear  form  associated  with  the  quadratic  form  S. 
Then  the  extremal  uq  of  the  quadratic  form  Q(u,  u,  A)  satisfies  the  equation 

Q(uo,  v,  A)  =  0,  Vu  G  Dq. 

Following  this  scheme,  from  Sobolev  embedding  theorems  we  obtain  the  fol¬ 
lowing  estimate  for  the  quadratic  form  Q(u,  u,  A): 

Q(u,  u,  A)|  <  C  ||  u  \\2Hl ,  C  =  const, 

which,  in  turn,  implies  the  following  representation  for  the  sesquilinear  form 
Q(u,  v,  A)  associated  with  Q(u,  u ,  A)  for  all  u,  v  G  Dq: 

Q(u,v,X)  =  (L(X)u,v)Hi,  (6) 

where  L( A)  :  H i(H)  — >  H i(H)  is  some  bounded  operator.  Let  us  find  the  form  of 
this  operator. 

Theorem  1.  The  operator  L( A)  :  -Hi(fl)  — >  H i(H),  associated  with  the  quadratic 
form  Q{u,  u,  A)  ,  has  the  form 

n 

L(X)u  =  u  —  (A  +  1)(Aq  +  I)~1u  +  —  Xmj)go(M ;  Af;  —  l)u(Mj), 

3= 1 

where  go(M ;  5;  —1)  is  the  Green’s  function  of  the  boundary  value  problem 

(A l  +  I)u  =  f(u ),  M  €  n, 

l(u)\on  =  0. 
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Proof.  Let  us  consider  the  sesquilinear  forms  bounded  on  H\  and  their  rep¬ 
resentations 

J  A0uA0vdxdy  =  (. Aiu,v)Hl , 
n 

J  uvdxdy  =  (. A2u ,  v)hx  , 
n 

u(M*)v(M *)  =  (A3(Msf)u,v)Hl, 

and  find  the  corresponding  operator  Ak,  k  =  1,2,3  for  each  form.  Since  all 
operators  A ^  :  H\  — »  H\  are  bounded,  then  it  suffices  to  consider  them  on 
functions  of  the  space  Cg°(J2). 

Note  that 

J  A0uA0vdxdy  =  (A0u,  A0v)L2^  =  (A lu,v)L^ay 
n 

On  the  other  hand,  for  the  same  functions  u,v  G  CQ°(f2)  we  have 
(Ani,v)Hl  =  (A0(Aiu),  A0v)L2(q)  +  (^iu,u)i2(n)  =  ((A§  + J)(Aiu),u)i2(n}, 
which  implies  the  validity  of  the  following  identity: 


((ao  +  /)(^i«)^)l2(Q)  =  (Aq  «,u)i2(n). 

Thus,  the  required  operator  has  the  representation 

Aw  =  (Ajj  +  (7) 

In  view  of  the  boundedness  of  the  operator  A\  :  Hi  — >  H i,  its  domain  is  the  space 
Hi. 

Reasoning  similarly,  from  the  identity 

J  uvdxdy  =  (A2u,v)Hl(n), 
n 

we  obtain  the  representation 


A2u  =  (Aq  +  I )  1u.  (8) 

Consider  the  third  identity  for  an  arbitrary  point  M*  =  M(x*,r/*)  G  12.  The 
following  lemma  is  true. 
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Lemma  1.  For  VM*  £  17,  it  £  H i(f2) 

A3(M*)u  =  go(M,  M*,  -l)u(M*).  (9) 

Proof.  It  is  well  known  that 

9q{M,  M*,  -1)  =  C\p2(M ,  M*)  In  p(M,  M*)  +  £(M),  C\  =  const, 
where  f(M)  is  such  that 

(A§  +  J)£(M)  =  0. 

Consequently,  (O). 

Front  the  embedding  theorem  for  Sobolev  spaces  we  obtain  the  estimate 

KADI  <||  u  ||c(n)<  b  II  u  ||/rl(n), 

where  6  >  0  is  a  constant  which  does  not  depend  on  the  choice  of  the  function 
it  £  H±  (fl) .  Then  A3(M*)  is  a  bounded  one-dimensional  operator.  It  remains  to 
show  that  for  all  u,v  £  Co°(fl)  the  identity 

u(M*)u(M*)  =  (. A3{M*)u,v)h i 

holds.  We  have 

{A3(M*)u,v)Hl  =  (A0(A3(M*)u),  A0u)L2(n)  + 

+(A3(M*)u,  u)x,2(n)  =  (A3(AQu,  (Aq  +  I)v)L2^y 
It  is  obvious  that  the  last  scalar  product  can  be  represented  in  the  form 

it(Ai*)  J  g0(M*;£,ir,-l)(Al  +  l)v(£,ri)d(;dri  = 
n 

=  u(  M*)(Aq  +  /)_1(Aq  +  /)u|m=m,  =  it(  AQv(AQ, 

which  proves  the  validity  of  the  representation  for  the  operator  A3(M*). 

The  lemma  is  proved.  ◄ 

From  formulas  (7)  -  (9)  we  obtain  the  following  representation  for  L( A): 

Q(u,  V,  A)  =  ( L(X)u ,  V )  =  (Am,  v)m  ~  X(A2u,  v)Hl  + 

n 

+  -  Xmj)(A3(Mj)u,v)Hl,  Vu,v  £  D q, 

3= 1 
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from  which  the  assertion  of  the  theorem  follows.  ◄ 

Note  that  the  constructed  operator  L( A)  :  Hi  — *  H\  has  an  obvious  connection 
with  the  biharmonic  differential  operator  pencil  with  5-like  coefficients. 

Indeed,  we  have  the  identity 

((Ao  +  v ) La(n)  =  (A 0u,  A0u)i2(Q)  -  A (u,  v)Lz(n)+ 

n 

+  ^2(kj  —  \rrij)8j(M  —  Mj)u(Mj)v(Mj),  Vu,v  €  Dq.  (10) 

3= 1 

The  closure  of  the  right-hand  side  of  (10)  defines  a  biharmonic  operator  with 
5-like  potential 


Lsu  =  (Aq  —  A  I)u  +  —  \rrij)5j(M  —  Mj)u(M) 

3=1 

with  domain  Dls  C  Dq. 

Questions  of  the  correct  determination  of  similar  operators  with  singular  co¬ 
efficients  were  studied,  for  example,  in  papers  [4?  ]. 

3.  Direct  spectral  problem 

Now  we  proceed  to  the  study  of  the  spectral  problem  and  its  properties: 

L(X)u  =  0,  u  €  H\,  A  6  C.  (11) 

The  solutions  of  the  spectral  problem  (11)  are  natural  frequencies  of  oscillations 
and  forms  of  natural  oscillations.  First,  we  show  that  the  spectrum  L( A)  is 
discrete  with  a  condensation  point  at  infinity. 

Denote 

Dq(o)u  =  (Aq  +  a/)-1,  u  G  Hi, 

Pj(a)u  =  go(M,  Mj, —a)u(Mj),  j  =  l,n,  g(a)  =  X  +  a. 

Let  us  note  that  if  —a  is  not  an  eigenvalue  of  the  operator  Aq,  then  the 
operator  (Aq  +  al)~1( Aq  +  I)  is  bounded  and  boundedly  invertible  in  the  space 

th¬ 
in  this  case  it  is  evident  that  the  equations  L(X)u  =  0  and  (Aq  -\-aI)~1( Aq  + 
I)L(X)u  =  0  are  equivalent.  Let  us  transform  the  last  equation: 

n 

D0(a)  +  ^2  kjPjipi) 

3=1 


u(M)  =  n(a) 


u— 
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-£(*=  j  +  rrij)Pj{a)u.  (12) 

3= 1 

In  this  equation,  the  operator 


D\(a)  =  Dq{q)  +  mjPj{a)  :  H\  -»  H\ 

3= 1 


is  compact,  and  the  operator 


£>2 («)  =  +  mj)pj(a)  =  ->• 

l=i 

is  bounded  and  finite-dimensional. 

Note  that  since  ||  D-2(a)  ||— >  0,  for  a  ->  +oo  and  a£l,  then  without  loss  of 
generality  we  can  assume  that  the  operator  (J  +  D2(a))~1  H\  H\  is  bounded 
for  any  set  of  fixed  parameters  kj  €  C ,nrij  G  C,  j  =  1  ,n.  Consequently,  equation 
(12)  can  be  rewritten  in  the  form 

u(M )  =  fj,(a)  (I  +  D2(a))~1  D\(a)u.  (13) 

Let  us  show  that  the  spectrum  of  the  operator  (/  +  D2{a))~l  D\(a)  consists 
of  an  infinite  series  of  eigenvalues  <7^  =  1  / fik(a)  with  a  single  condensation  point 
at  the  zero  point. 

The  operator  ( I  +  D2(a))~1  D\{a)  :  H\  H\  is  compact,  hence  its  spectrum 
consists  only  of  eigenvalues  of  finite  multiplicity  except  possibly  a  zero  point. 

It  is  easy  to  show  that  the  dimension  of  the  kernel  of  the  operator  (I  + 
T>2(a)) ~1Di(a)  can  not  be  greater  than  n  +  1,  and  its  image  is  an  infinite¬ 
dimensional  subspace  of  ITf  (fi). 

Thus,  the  following  theorem  is  true. 

Theorem  2.  The  spectrum  of  the  operator  pencil 

n 

L(X)u  =  u  —  (A  +  1)  (/  +  Aq))  1  u  +  ^2(kj  —  Xmj)go(M,  Mj ,  —  l)u(Mj) 

3= 1 

for  an  arbitrary  set  of  parameters  k  =  (ki,..,kn)  C  Cn,rti  =  (mi,..,m„)  C  Cn 
is  discrete  with  a  single  condensation  point  at  infinity. 
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Let  A  be  a  regular  point  of  the  unperturbed  operator  Lq(X)  :  H\  — >  H\,  where 

L0(A)u  =  it  —  (A  +  1)(7  +  Aq)_1u.  (14) 

Then  the  operator  (Ag  —  AI)_1(I  +  Ag)  is  boundedly  invertible  in  H\  and  the 
spectral  problems  L(X)u  =  0  and  L{X)u  =  (Ag  —  XI)~l(I  +  Ag)L(A)n  =  0  are 
equivalent.  It  is  clear  that  the  equation  L(X)u  =  0  can  be  written  in  the  form 

n 

L{X)u  =  u(M,  A)  —  ^~~](A nij  —  kj)go(M ,  Mj,  X)u(Mj ,  A)  =  0.  (15) 

3= 1 

Considering  the  equation  (15),  let  us  show  that  only  at  the  points  M  = 
Mj ,  j  =  l,n,  it  is  equivalent  to  a  finite-dimensional  system.  For  this  purpose,  let 
us  write  the  equation  (15)  for  a  series  of  points  M  =  Mj ,  j  =  1  ,n: 

u(Mi,  A)  =  Y!j=i(^mj  ~  kj)g0(Mi,  Mj,  X)u(Mj,  A) 
u(M2,  A)  =  Z)"=1(A mj  -  kj)g0(M2,Mj,X)u(Mj,X) 

,  u(Mn,  A)  =  YIj=i^mj  ~  kj)g0(Mn ,  Mj,  X )u(Mj,  A) 

Let 


(  5o (Ml,  Mi,  A)  ...  go(M\ ,  ALn,  A)  \ 

/o  go(M2,  M\ ,  A)  ...  go(M2,  Mn,  A) 

\  go(Mn,  Mi,  A)  ...  go(Mn,Mn,X)  J 

A  =  diag{k\,  k2, ...,  /cn},  B  =  diag{m\,m2 , ...,  mn}  are  diagonal  matrices  and 
u(A)  =  (u(Mi,  A),  u(M2,  A),  ...,u(Mn,  X))T .  Then  the  series  of  obtained  equalities 
can  be  written  in  the  following  matrix  form: 

(/  +  G0(A)(A-A5))u(A)  =  0.  (16) 

For  the  solvability  of  system  (16)  it  is  necessary  and  sufficient  that 

/(A)  =  det(/  +  G0(A)(A-AB))  =  0.  (17) 

Let  us  note  that  the  meromorphic  function  /(A)  with  poles  at  the  points  A  G  cx(Lo) 
is  the  spectrum  of  the  (discrete)  operator  Lq. 

Let  X*  be  a  zero  of  the  function  /(A).  Then  it  follows  from  (15)  that  the 
eigenfunctions  of  the  operator  L( X)  admit  a  representation 

n 

u(M,  A*)  =  £(A *mj  -  kj)go(M,  Mj,  X*)u(Mj,  A*). 


(18) 
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Thus,  the  following  statement  is  true. 

Theorem  3.  Each  eigenvalue  of  the  pencil  L{ A)  for  any  set  of  parameters  k  = 
(ki,..,kn)  C  C  n,rfr  =  (mi, ..,  mn)  C  Cn  is  either  a  zero  of  f(  A),  or  is  an  eigen¬ 
value  of  the  pencil  Lq(\). 

4.  The  inverse  spectral  problem  and  its  properties 

Consider  the  inverse  spectral  problem  for  the  operator  pencil  L( A),  which 
linearly  depends  on  the  parameters  k  =  (k\, .. kn ),  m  =  (mi,  ..,mn ): 

L(\)u  =  L( A,  k,  m)u  = 

n 

=  u  -  (A  +  1)  (I  +  Aq))  1  u  +  y ^(kj  -  \mj)g0(M,  Mj,  -1  )u(Mj). 

3= 1 

The  eigenvalues  of  the  pencil  L( A,  k,  m)  also  depend  on  the  parameters  k  = 
(ki, .. kn ),  th  =  (mi,  ..,mn),  and  this  dependence  is  analytic,  with  the  exception 
of  singular  points  that  form  a  set  of  measure  zero. 

Let  us  formulate  a  multiparameter  inverse  spectral  problem  (MPISP)  for  the 
operator  pencil  L(\,k,m). 

MPISP  It  is  required  to  choose  such  values  mj  of  point  masses  and  stiffness 
coefficients  of  springs  kj,  that  2 n  of  eigenvalues  of  the  operator  pencil  L( A,  k,  m) 
are  equal  to  the  given  numbers  A^,  A?;, ...,  A* ,  A*+1, ...,  A|n . 

In  this  case  A*  =  (A|,  A?;, ...,  A*  ,  A*+1, ...,  A^)  G  C2n  will  be  called  the  spectral 
data  of  MPISP.  By  the  solution  of  MPISP  we  will  mean  a  set  of  parameters 
k,  m  €  Cn,  such  that  dim.  ( ' ker(L(\ *,  k,  m))^  >  0  for  Vj  =  1,  2 n. 

Further  we’ll  show  that  MPISP  for  the  operator  pencil  L( A,  k,  m)  can  be  rep¬ 
resented  as  MPISP  for  a  finite-dimensional  operator.  Let  us  rename  the  matrices 
and  parameters  for  convenience. 


mi  =  pi,m2  =  p2-..,mn  = 

=  Pn, 

k{  —  Pn- 1-1)  •••)  kn  —  pi, 

1  = 

2  n 

(  0 

...  go(M\ ,  Mk,  A)  ... 

0  ^ 

B0(X)  =  I,  Bk(  A)  = 

0 

...  go(M2,Mk,\)  ... 

0 

V  o 

...  g0(Mn,Mk,X )  ... 

0  ) 

Bn+k( A)  =  A  •  Bk(X)  k  =  1,  n. 
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Then  the  system  (16)  can  be  written  in  the  form 

B(p,  X)u(X)  =  (Bo(X)+p1B1(X)  +  ...+plBl(X))u(X)  =  0,  (19) 

where  p  =  (pi,p2,  Pi )  G  Cl. 

Note  that  B(p,  A)  :  C  — >  En  for  each  fixed  vector  p  is  a  meromorphic  function 
in  A  G  C  with  poles  at  the  points  of  the  spectrum  of  the  unperturbed  pencil  Lq(X). 
Thereby,  we  will  further  require  that  the  spectral  data  A*  =  (A*,  A|, A)1)  G  of 
MPISP  do  not  intersect  with  the  spectrum  of  an  unperturbed  pencil  of  operators 

L0(X). 

Let  us  reformulate  the  initial  MPISP  for  an  infinite-dimensional  pencil  L( A,  k,  m) 
H i  — >  H i  equivalent  to  the  MPISP  for  a  finite-dimensional  matrix  pencil  B(p ,  A)  : 
En  -a  En. 

MPISPf  It  is  required  to  find  possible  values  of  the  vector  p  from  the  space 
Cl  ,  for  which  the  given  numbers  Ai,  A2, A;  are  the  eigenvalues  of  the  operator 
B(p,X),  i.e.  det(B(j>,  A))  =  0. 

Thus,  under  the  assumption  that  X*k  is  not  an  eigenvalue  of  the  pencil  Lo(A), 
the  inverse  spectral  problem  for  the  operator  pencil  is  reduced  to  a  system  of 
l  =  2 n  polynomial  equations  of  l  =  2n  variables  pi,  ..,pi 

det(B(p,  A*))  =  0,  j  =  1,  l. 

It  can  be  argued  that  the  solution  of  the  inverse  spectral  problem  for  L( A,  k,  m) 
is  the  union  of  algebraic  manifolds  in  C;  and  a  finite  set  of  points  in  C;  . 

The  existence  and  isolation  problems  of  the  solutions  of  MPISP  for  oscillations 
of  elastic  plates  and  shells  with  dotted  fastenings  will  be  discussed  in  detail  in  a 
separate  paper. 

In  conclusion,  note  that  the  mathematical  model  proposed  in  this  paper  for 
describing  an  elastic  plate  with  point  fixings  can  also  be  used  to  study  the  elastic 
oscillations  of  beams,  shells,  and  bulk  bodies  with  additional  point  fixings. 
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